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I Abstract 

In this paper we first give a short overview of the MACRO detector which was operational at the Gran 
Sasso Laboratory from 1989 till the end of 2000. Then we present and discuss the results on atmospheric 
muon neutrino oscillations, concerning medium (~ 4 GeV) and high (~ 50 GeV) energy data. Using 
the Multiple Coulomb Scattering of muons inside the lower part of the detector, estimates of neutrino 
I energies were made for the high energy sample. The data on angular distributions, absolute flux and L/E^, 
' distributions favour i/^ < — > 1/^- oscillations with maximal mixing and Ato^ = 2.5 • 10^"^ cV^. 



X ' 1 Introduction 

V ' 

Qh! macro was a large multipurpose underground detector designed primarely to search for 
^ I rare particles in the penetrating cosmic radiation. Though it was optimized to search for 
the supermassive magnetic monopoles predicted by Grand Unified Theories (GUT) of the 
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electroweak and strong interactions, the MACRO research program included a wide range 
of topics in areas of astrophysics, nuclear and particle physics and cosmic ray physics. In 
I particular it planned to study atmospheric neutrinos via the detection of upgoing muons 
produced in charge current (CC) interactions of in the rock below the detector and inside 
its lower part. 

A high energy primary cosmic ray, proton or nucleus, interacts in the upper atmosphere 
producing a large number of pions and kaons, which decay yielding muons and muon 
neutrinos; also the muons decay yielding muon and electron neutrinos. Fig: 0. The ratio of 
the numbers of muon to electron neutrinos is about 2, and N^/Ny ~ 1. These neutrinos 
are produced in a spherical surface at about 20 km above ground and they proceed at high 
speed towards the Earth. An underground detector is "illuminated" by a flux of neutrinos 
from all directions. 

At low energies, ^ 1 GeV, the numbers of atmospheric neutrinos predicted by different 
authors differ by about 20 -i- 30% |l[]; at higher energies, Ey > 10 GeV, the predictions are 
more reliable, with an estimated systematic uncertainty of about 15% 0, almost one half 

^ For the full list of the Collaboration see the last paper of ref . |b|l 
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Figure 1: (a) Illustration of the production, travel and interactions of atmospheric muon neutrinos; 
(b) interaction of a primary cosmic ray, production of pions (and kaons) and their decays leading to the 
atmospheric Ve, t'/j- 



of the uncertainty at low energies. However the predicted relative rates of to z/g and the 
shapes of the zenith angle distributions are affected by considerably lower systematic errors. 
Other sources of systematic uncertainties arise from the knowledge of the neutrino-nucleon 
cross sections and from the propagation of muons and electrons in different materials. 

Several large underground detectors studied atmospheric neutrinos, mainly via charged 
current interactions, where a z/^ gives rise to a and thus to a track, and a Ug yields an 
e~ and thus to an electromagnetic shower. 

The early water Cherenkov detectors 1MB and Kamiokande reported anomalies in 
the double ratio of muon to electron neutrinos integrated over all zenith angles, while the 
tracking calorimeters NUSEX |^ and Frejus @], and the Baksan scintillator detector did 
not find any deviation, see Fig: ^. 

Later the Soudan 2 p|, MACRO and SuperKamiokande |]10[ detectors reported a 
deficit in the z/^ flux with respect to the MC predictions and a distortion of the angular 
distributions, while the z/g flux and angular distributions agree with MC. These features 
may be explained in terms of z/^ < — > Vr oscillations. 

Atmospheric neutrinos are well suited for the study of neutrino oscillations, since they 
have energies from a fraction of a GeV up to more than 100 GeV and they may travel 
distances L from few tens of km up to 13000 km; thus the LjEy values range from ~ 1 
km/GeV to more than 10*^ km/GeV. 

In the following, we shall briefly recall neutrino oscillation formulae and the main features 
of the MACRO detector, then we shall discuss the main results of MACRO, including the 
ratio of vertical to horizontal muons which favour z/^ < — > oscillations and the estimates 
of neutrino energies by multiple coulomb scattering of muons, which allows the study of the 
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Figure 2: Double ratios i?' measured by several atmospheric neutrino experiments. 

L/Ey distribution. We shall conclude with a short discussion and by recalling briefly other 
neutrino studies by MACRO. 



2 Neutrino oscillations 

If neutrinos have non-zero masses, one has to consider the neutrino weak flavour eigenstates 
Vei i^ti, and the mass eigenstates ui, 1/2, 1^3. The weak flavour eigenstates z/; are linear 
combinations of the mass eigenstates through the elements of the mixing matrix Uim'- 
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^1= ^Irn l^m (1) 

m=l 

If the mixing angles are small, one would have i/g ~ '^i? ~ ^^2? '^r ~ ^'a- If the mixing 
angles are large, the flavour eigenstates are well separated from those of mass. 

In the simple case of only two flavour eigenstate neutrinos (i/^, Ur) which oscillate with 
two neutrino mass eigenstates (z/2, z/3) one has 

= z/2 cos 023 + V3 sin 623 
Vt = -V2 sin Q23 + cos Q23 

where Q23 is the mixing angle. In this case one may easily compute the following expression 
for the survival probability of a beam (disappearance experiment): 

P{v, - = 1 - sin^ 2^23 sin^ (^^i) = 1 - sin^ 2^23 sin^ |^ l-27Am^-^ ^ 

where Am^ = — vp?^^ and L is the distance travelled by the neutrino from production to 
the detection point. The probability for the initial z/^ to oscillate into a is (appearance 
experiment): 

P(/^M ^i^r) = l- P{y>. ^ /^m) = sin' 2^23 sin^ ( — ^ j (4) 

In a disappearance experiment one measures only P{v^ — > z^^) as a function of Ey and 
L, or of L/Ey. Only disappearance experiments have been performed until now. For 



3 



(a) 



(b) 



Figure 3: (a) Vertical section of the MACRO detector. Event topologies induced by interactions in or 
around the detector: lU = semicontained Internal Upgoing /i; ID — Internal Downgoing fi; UGS — Upgoing 
Stopping /J,; Upthroughgoing = upward throughgoing /i. The black circles indicate the streamer tube hits 
and the black boxes the scintillator hits, (b) MC simulated distributions of the parent neutrino energies, 
assuming no oscillations, giving rise to the different topologies of muon events detectable by MACRO. The 
distributions are normalized to one year of data taking. 

atmospheric neutrinos the parameters 623 and Am^ may be determined from the variation 
of P{iy^ — > i^fj.) as a function of the zenith angle B, or from the variation in L/E^. In future 
appearance experiments, one would try to observe a neutrino, for example Ur, not produced 
in 71, K decays. 

3 The MACRO experiment 

The MACRO detector was located in Hall B of the Gran Sasso Laboratory, under a minimum 
rock overburden of 3150 hg/cm^ which reduces the atmospheric muon flux by a factor 
~ 5 ■ 10^. It was a large rectangular box, 76.6mx 12mx9.3m, divided longitudinally in six 
supermodules and vertically in a lower part (4.8 m high) and an upper part (4.5 m high) 
(Fig: ^). It had three types of detectors which gave redundance of informations: liquid 
scintillation counters, limited streamer tubes and nuclear track detectors. This last detector 
was used only for new particle searches. 

For muon and for neutrino physics and astrophysics studies, the streamer tubes were 
used for muon tracking and the liquid scintillation counters for fast timing. The lower part 
of the detector was filled with trays of crushed rock absorbers alternating with streamer 
tube planes; the upper part was open and contained the electronics. 

There were 10 horizontal planes of streamer tubes in the bottom half of the detector, and 

4 planes at the top, all with wire and 27° stereo strip readouts. Six vertical planes of streamer 
tubes and one layer of scintillators covered each side of the detector. The scintillator system 
consisted of three layers of horizontal counters, and of the mentioned vertical layer along 
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Figure 4: distributions for throughgoing (downgoing and upgoing) muons for the runs with the complete 
detector. The shaded areas concern muons crossing 3 scintillation counters, while the open histogram 
concerns muons crossing 2 or 3 counters. Two vertical dotted lines delimit the range —1.25 < < —0.75 
of upthroughgoing muons. There are ^ 3.5 • 10^ downgoing muons with 1//? ~ 1 and 782 (before subtraction 
of background) upgoing muons with 1//3 ^ —1. 

the sides of the detector. One thus had a closed box structure with openings only in the 
upper part of the front and back ends. 

The time and space resolutions for muons in a scintillation counter were about 500 ps and 
11 cm, respectively, while in the streamer tube system they were about 100 ns and 1 cm, 
respectively. The combination of the informations from the streamer tubes and from the 
scintillators allowed tracking with a precision of 1 cm over path lengths of several meters, 
and timing with a precision of 500 ps. The detector provided a total acceptance S'fi ~ 10000 
m^ sr for an isotropic flux of particles. 

Fig: l^a is a vertical section of the detector; it shows a general view of the detector and 
gives also a sketch of the different topologies of detected neutrino-induced muon events 
used to study neutrino oscillations: Upthroughgoing muons. Internal Upgoing muons (lU), 
Upgoing Stopping muons (UGS) and Internal Downgoing muons (ID). Fig: ^3 shows the 
parent neutrino energy distributions obtained from Monte Carlo calculations for the event 
topologies detectable in MACRO. 

The Upthroughgoing muons (with E^> 1 GeV) come from interactions in the rock below 
the detector of muon neutrinos with an average energy {Ey) ~ 50 GeV. The tracking is 
performed with streamer tubes hits; the time information, provided by scintillation counters, 
allows the determination of the direction (versus) by the time-of-flight (T.o.F.) method. If 
the atmospheric neutrino anomalies are due to neutrino oscillations, one expects a muon 
flux reduction depending on the zenith angle. 

The semicontained upgoing muons (lU) come from interactions inside the lower 
apparatus. Since two scintillation counters are intercepted, the T.o.F. method is applied 
to identify the upward going muons. The average parent neutrino energy for these events 
is ~ 4 GeV. If the atmospheric neutrino anomalies are the results of oscillations with 
maximal mixing and 10^^ < Aiv? < 10~^ eV^, one expects a reduction of about a factor of 
two in the flux of lU events, without any distortion in the shape of the angular distribution. 

The up stopping muons (UGS) are due to interactions in the rock below the detector 
yielding upgoing muon tracks stopping in the detector; the semicontained downgoing muons 
(ID) are due to induced downgoing tracks with vertex in the lower MACRO. The 
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Figure 5: (a) Zenith distribution of the upthroughgoing muons (809 events, background subtracted). The 
data (black points) have statistical and systematic errors added in quadrature. The shaded region indicates 
the theoretical scale error band of ±17% (see text). The solid line is the fit to an oscillated flux with maximal 
mixing and Am^ = 2.5 • 10^"^ eV^. (b) Zenith distribution of the ratio of measured upgoing muons over the 
expected ones with no oscillations (black points with statistical and systematic errors in quadrature). The 
solid line is the same ratio computed for oscillations; the dotted line corresponds to no oscillations. 

events are found by means of topological criteria; the lack of time information prevents 
to distinguish between the two subsamples. An almost equal number of UGS and ID events 
is expected. In case of oscillations, a similar reduction in the flux of the up stopping events 
and of the semicontained upgoing muons is expected; no reduction is instead expected for 
the semicontained downgoing events (which come from neutrinos which travelled ~ 20 km). 

4 High energy events 

The data were collected during the running period from March 1989 to April 1994 with 
the detector under construction and during the runs with the complete detector from 1994 
until December 2000 (livetime 5.52 yrs). Since the total livetime normalized to the full 
configuration is 6.16 yrs, the statistics is largely dominated by the full detector run. The 
analysis of a data sample of more than 40 million atmospheric downgoing muons achieved 
a rejection factor of ~ 10"'' which includes background caused by showering events and 
radioactivity in coincidence with muons. 

We studied many sources of possible background: (i) low energy upgoing particles (mainly 
pions) induced by undetected downgoing muons in the rock surrounding the detector |I2| . 
This source of background is reduced to 1% requiring that more than 200 gr/cm^ of material 
be crossed in the lower part of MACRO; (ii) we excluded a region in the azimuthal angle 
from —30° to 120° for nearly horizontal upgoing muons due to insufficient rock overburden; 
(iii) for muons crossing 3 scintillation counters a linear fit of the times as a function of the 
path length is performed and a cut is applied on the x^; (iv) further minor cuts are applied 
to events crossing 2 counters. 
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Figure 6: 90% C.L. allowed region contours for z^^ < — > Vr oscillations obtained by the SuperKamiokande, 
MACRO and Soudan 2 experiments. 



Muons traveling downwards have ~ +1, upwardgoing muons have ~ —1. The 
distribution for the sample collected with the full detector is shown in Fig: ^ We 
selected upwardgoing muons requiring —1.25 < < —0.75; we found 863 events, of which 
809 events remain after subtracting the background. 

In the upgoing muon simulation the neutrino flux computed by the Bartol group |T^ was 
used. The cross sections for the neutrino interactions was calculated using the deep inelastic 
parton distributions The muon propagation to the detector was computed using the 



energy loss calculation in standard rock [|15|] . 

The systematic uncertainty arising from the neutrino flux, cross section and muon 
propagation on the expected flux of upthroughgoing muons was estimated to be ~17%. 
This systematic uncertainty is mainly a scale error that doesn't affect the shape of the 
angular distribution. The same cuts applied to the data were used for the simulated events: 
they selected 1122 MC events, assuming no oscillations. 

Fig: ^ shows the zenith angle (G) distribution of the measured flux of upthroughgoing 
muons of energy ii^^ > 1 GeV (black points); the MC expectation for no oscillations is 
indicated by the dashed line with the shaded scale error band. A deficit of events in the 
region around the vertical can be noticed. The ratio of the observed number of events 
to the expectation without oscillations in —1 < cos© < (B is the zenith angle) is 
0.721 ± 0.026,t,t ± 0.043,,, ± 0.123^;,. 

We tested the shape of the angular distribution with the hypothesis of no oscillations 
normalizing the prediction to the data. The jD.o.F. corresponds to a probability of 0.2%. 

In the hypothesis of only z/^ < > oscillations, the minimum x^/D.o.F. in the physical 

region is 9.7/9 (P = 37%) for maximal mixing and Am^ = 2.5 ■ 10~^ eV^. An independent 
test was made on the number of the events. 

Combining the probabilities from the two independent tests on the zenith angle shape of 
the flux and on the total number of events, the maximal probability is 66% for maximum 
mixing and Am^ ~ 2.4 ■ 10"'^ eV^. The result of the fit is the solid line in Fig: |^. 

For each bin of Fig: |^a we computed the ratio R = (measured muon flux) / (expected 
muon flux in the non oscillation hypothesis). The result is shown in Fig: (black points). 
The solid line is the ratio R for maximal mixing and Am^ = 2.5 ■ 10~^ eV^; the dashed line 
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Figure 7: Ratio of events with —1 < cos 9 < —0.7 to events with —0.4 < cos 9 < as a function of Am^ 
for maximal mixing. The black point with error bar is the measured value, the solid line is the prediction 
for < — > oscillations, the dotted line is the prediction for i'^ < — > t^sterUe oscillations. 



is the same ratio expected for no oscillations. 

The MACRO upthroughgoing muons are very sensitive to oscillations because of the 
energy range covered. The MACRO 90% C.L. allowed region in the sin^ 29 — Am^ plane, 



computed using the Feldman-Cousins method |16], is smaller than the SuperKamiokande 
(SK) one for upthroughgoing muons due to the different energy thresholds (~ 1 GeV for 
MACRO and ~ 6 GeV for SK) and for the median energies (~ 50 GeV for MACRO, 70^80 
GeV for SK). 

The MACRO 90% C.L. allowed region for < — > v^. oscillations is compared in Fig: |^ 
with those obtained by the SK and Soudan 2 experiments. 



5 Matter effects, < > against i/^ ^ > v sterile 

Matter effects due to the difference between the weak interaction effective potential for 
muon neutrinos with respect to sterile neutrinos (which have null potential) would produce 
a different total number and a different zenith distribution of upgoing muons. 

In Fig: ^the measured ratio between the events with —1 < cos© < —0.7 and the events 
with —0.4 < cosO < is shown by the black point. In this ratio most of the theoretical 
uncertainties on neutrino flux and cross sections cancel. The remaining theoretical error, 
estimated at < 5%, has a 3% contribution due to uncertainties on the kaon/pion fraction 
produced in atmospheric showers, a 2% contribution from uncertainties in the cross sections 
of almost vertical and almost horizontal events and a 2.5% contribution from the seasonal 
variations of the ratio. The systematic experimental error on the ratio, due to analysis 
cuts and detector efficiencies, is 4.6%. Combining the experimental and theoretical errors 
in quadrature, a global extimate of 7% is obtained. MACRO measured 305 events with 
-1 < cosO < -0.7 and 206 with -0.4 < cos0 < 0; the ratio is = 1.48±0.13stat±0.10,y,. 
For Am^ = 2.5 ■ 10"'^ eV^ and maximal mixing, the minimum expected value of the 
ratio for t'^ < — > Vr oscillations is Rr = 1.72; for z/^ < — > v sterile oscillations one expects 
Rsteriie = 2.16. The maximum probabilities Pbest to find a value of Rr and of Rsteriie smaller 
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Figure 8: Measured zenith distributions (a) for the lU events and (b) for the ID+UGD events. The black 
points are the data, the shaded regions correspond to MC predictions assuming no oscillations. The full 
line is the expectation for < — > Vt oscillations with maximal mixing and Am^ = 2.5 x 10"'^ eV^. 



than Rexpected are 9.4% and 0.06% respectively. Hence the ratio of the maximum probabihties 
is Pbestr / Pbeststertie ~ 1^'''' ^^^^ ^ ^^ ^ — ^ ^sterile osciUations (with any mixing) are excluded 
at 99% CL. compared to the < — > Uj. channel with maximal mixing and Am^ = 2.5- 10~^ 



6 Low energy events 

The data concern only the running period with the detector in the full configuration from 
April 1994 to December 2000. During this period more than 40 million downgoing muons 
were collected. Because of the difference between the topologies of the low energy events, 
two separate analyses were performed. 

The lU sample corresponds to an effective livetime of 5.8 yrs. The basic request is the 
presence of a streamer tube track reconstructed in space matching at least two hits in two 
different scintillators in the upper part of the apparatus. The measured muon velocity (3c was 
evaluated with the sign convention that upgoing (downgoing) muons have ~ — 1 (~ +1). 
After background subtraction, we had 154 upgoing partially contained events. 

The identification of ID+UGS events was based on topological criteria. The candidates 
had a track starting (ending) in the lower apparatus and crossing the bottom detector face. 
The track has also to be located or oriented in such a way that it could not have entered 
(exited) undetected through insensitive regions of the apparatus. For this analysis the 
effective livetime was 5.6 yrs. After background subtraction, we had 262 ID+UGS events. 
The MC simulations for the low energy data use the Bartol neutrino flux [r^| and the 



low energy neutrino cross sections |]T^. We estimated a total theoretical scale uncertainty 
on the predicted number of muons of the order of 25% (it is very probably overestimated). 

With the full MC simulation, the prediction for lU events was 285 ± 28syst ± 7^th, while 
the observed number of events was 154 ± 12stat- The ratio was Riu = (Data/MC)/;/ = 
0.54 ± O.Oistat ± 0.05syst ± 0.13th- For this sample, the average value of log^^Q L/Ei, is 3.2, in 
good agreement with the hypothesis of neutrino oscillations. 

The prediction for ID+UGS events was 375 ± 37syst ± 94t/i, while the observed number 
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Upthroughgoing 



lU 
ID+UGS 



Events 
809 
154 
262 



MCno osc 
1122 ± 191 
285 ± 28,^,, ± 71th 
375 ± 37,ys ± Mth 



R = (Data/MCno osc) 
0.721 ± 0.026,t,t ± 0.043,^, ± 0.1234^ 
0.54±0.04,t,t±0.05,j,,±0.13tft 
0.70 ±0.04,t,t± 0.07,2,, ±0.17t,, 



Table 1: Summary of the MACRO muon events in — 1 < cos0 < after background subtraction. For 
each topology the number of measured events, the MC prediction for no oscillations and the ratio R = 
(Data/MCno osc) are given. 



of events was 262 ± IGstat- The ratio was Rid+ugs = (D^'^^/^G)j£)+ugs = 0.70 ± O.O^stat ± 



The numbers of low energy events are compared with the predictions in Table: ^ the 
angular distributions are shown in Fig: |^. The data show a uniform deficit of the measured 
number of events over the whole angular distribution with respect to the predictions without 
oscillations. The two data sets are consistent with neutrino oscillations with the parameters 
found in the analysis of the high energy sample. Upgoing neutrinos which induce lU and 
UGS events, travelling thousands of kilometers, are reduced by 50%, while no reduction is 
expected for downgoing partially contained muons. As a rough prediction, we expect a rate 
reduced by 50% for lU and by 25% for ID+UGS. 

The measured data are in good agreement with the predicted angular distributions based 

on z/^ < y Ur oscillations with the parameters obtained from the upthroughgoing muon 

sample (full histogram in Fig: 

7 i^^ energy estimates by Multiple Coulomb Scattering of 
upthroughgoing muons 

Since MACRO was not equipped with a magnet, the only way to experimentally estimate 
the muon energy is through Multiple Coulomb Scattering (MCS) of muons in the absorbers 
in the lower part of MACRO, see Fig: The r.m.s. of the projected displacement of a 
relativistic muon with momentum p travelling for a distance X can be written as: 



where p (GeV/c) is the muon momentum and X/X^ is the amount of crossed material 
in units of radiation lengths. A muon crossing the whole apparatus on the vertical has 
aMcs — 10(cm)/£'(GeV). The muon energy estimate can be performed up to a saturation 
point, occurring when aMcs is comparable with the detector space resolution. 
Two analyses were performed. 

The first analysis was made studying the deflection of upthroughgoing muons in MACRO 
with the streamer tubes in digital mode. Using MC methods to estimate the muon energy 
from its scattering angle, the data were divided into 3 subsamples with different average 
energies, in 2 samples in zenith angle B and finally in 5 subsamples with different average 
values of L/E^. This method could reach a spatial resolution of ~ 1 cm. This analysis 
yielded a.n L/E,y distribution quite compatible with neutrino oscillations with the parameters 
found in Section 4. 



C^MCS - 



X 13.6 ■ IQ-^GeV 
V3 ppc 



^JX/X^ ■ (1 + 0.038 ln(X/X°)) 



(5) 
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Figure 9. R = (Data/MCno osc) as a function of \ogiQ{L/ E,y) obtained using the streamer tubes in "drift 
mode". The data are the points with error bars, the vertical extent of the shaded regions represents the 



uncertainties on the MC predictions for t/^ 



oscillations and for no oscillations. 



As the interesting energy region for atmospheric neutrino oscillations spans from ~ 1 GeV 
to some tens of GeV, it is important to improve the spatial resolution of the detector to push 
the saturation point as high as possible. For this purpose, a second analysis was performed 



using the streamer tubes in "drift mode" [|l8j using the TDC's included in the QTP system. 



originally designed for the search for magnetic monopoles. To check the electronics and 
the feasibility of the analysis, two "test beams" were performed at the CERN PS-T9 and 
SPS-X7 beams. For each triggered tube, the arrival time of the signal multiplied by the drift 
velocity, measured at the test beam, gives the radius of the drift circle. The upthroughgoing 
muon tracks were reconstructed as the best fit of the drift circles referring to the triggered 
streamer tubes. The space resolution achieved is ~ 3 mm, a factor of 3.5 better than in 
the previous analysis. For each muon, 7 MCS sensitive variables were given in input to a 
Neural Network (NN) which was previously trained with MC events of known input energy 
crossing the detector at different zenith angles. The NN output allows to separate the whole 
sample of upthroughgoing muons in 4 subsamples with average energies Efj^ of 12, 20, 50 
and 102 GeV, respectively. The comparison of the zenith angle distributions of the 4 energy 
subsamples with the predictions of no oscillations MC shows a strong disagreement at low 
energies (where there is a deficit of vertical events), while the agreement is restored at higher 
neutrino energies. The corresponding probabilities for the no oscillation hypothesis in 
these energy windows are 1.8%, 16.8%, 26.9% and 87.7%, respectively. 

The distribution of the ratios R = (Data/MCno osc) obtained by the second method is 
plotted in Fig: ^ as a function of logiQ^L / E^) . The points with error bars are the data; the 
vertical extent of the shaded areas represents the global uncertainties on the MC predictions 
(i) for < — > Uj. oscillations with maximal mixing and Am^ = 2.5 x 10~^ eV^ and (u) for 
no oscillations. 

The 90% C.L. allowed region from this analysis seems to be consistent with the MACRO 
region plotted in Fig: M. 
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8 Conclusions 



We presented and discussed the MACRO results on atmospheric neutrinos. The measured 
and expected numbers of events are summarized in Table: ^ The observed zenith angle 
distributions and the numbers of neutrino-induced muons (Fig: ^ and Fig: ^ disagree with 
the predictions of the no oscillation hypothesis. The muon angular distribution and flux are 

in agreement with the hypothesis of two flavour z/^ < > Vr oscillations, with maximal mixing 

and Am^ = 2.5 ■ 10~^ eV^. The hypothesis of < — > VsterUe oscillations is disfavoured at 
the 99% C.L. The MACRO 90% C.L. contours in the plane sin^ 26 - Am"^ are compared 
in Fig. ^ with those obtained by the Soudan 2 and SuperKamiokande experiments: they 
overlap closely. The data from the three experiments strongly favour z/^ < — > Uj. oscillations, 
though other possible interpretations have been suggested, see for example ref. [1T9| . 

Other MACRO results in the field of neutrino physics and astrophysics concern: (i) the 
search for cosmic point sources of high energy muon neutrinos: the flux upper limits for 
many sources, like Vela P, are at the level of $jy ~ 2 • 10~^ cm~^ s~^ ||20|; (ii) the search for 
dark matter WIMPs annihilating and yielding muon neutrinos from the earth core and from 
the Sun: the measured flux upper limits exclude ranges of parameters of supersymmetric 
models ; (iii) the search for low energy from stellar gravitational collapses: from 1989 
till 2000 no type II supernovae exploded in our galaxy [52 . 
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